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c^ ■ Abstract 

iu: 

r-| , We define a simplified version of Regge quantum gravity where the link lengths can take 



on only 2 possible values, both always compatible with the triangle inequalities. This is 
therefore equivalent to a model of Ising spins living on the links of a regular lattice with 



^ • somewhat complicated, yet local interactions. The measure corresponds to the natural 

sum over all 2*'*"'^'^ configurations, and numerical simulations can be efficiently imple- 
mented by means of look-up tables. In three dimensions we find a peak in the "curvature 
susceptibility" which grows with increasing system size. However, the value of the corre- 
sponding critical exponent as well as the behavior of the curvature at the transition differ 
from that found by Hamber and Williams for the Regge theory with continuously varying 
link lengths. 
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I. INTRODUCTION 



To date, two main formulations of lattice quantum gravity have been considered, the so 
called "Regge gravity" [|l], and "simplicial gravity" approaches 0, §, ||. While it could 
be argued that both formulations involve Regge calculus and simplexes, the distinguishing 
feature is that the former has a fixed incidence matrix and varying link lengths while the 
latter has a varying incidence matrix and fixed link lengths. 

Both formulations are technically and computationally demanding. For example, the 
Regge approach involves calculating areas and deficit angles involving general d-simplexes. 
In the simplicial approach these take on only a limited number of possible values, but 
the updating moves involve complicated interchanges of several simplexes at once. We 
introduce here a third lattice gravity approach which is structurally and computationally 
much simpler than either Regge or simplicial gravity, and may as a result be amenable to 
analytic attack in more than two dimensions. 

We call our formulation "Ising-link quantum gravity" . It is easy to define. The incidence 
matrix is fixed exactly as in the Regge approach. But the link lengths can only take on 

two values, 



li = l + bsi (1) 

with Si = ±1 and b a positive constant.^ i is a link label. In order that the triangle inequal- 
ity (or its higher- dimensional generalization - that the simplex volume is real and positive) 
is always satisfied, it is straightforward to show that we must take 6 < | in two dimensions, 
& < 3 — Vs ~ .17 in three dimensions, etc. (See Section II.) We restrict b to satisfy this 
inequality so that all 2^^ configurations are allowed. (A^i is the number of links.) This is 
quite different from either Regge or simplicial gravity where most potential updates either 
violate the triangle inequalities or violate the manifold property. Furthermore it provides 
us with a natural measure which gives all 2^^ configurations equal weight. It is clear that 
our model is completely equivalent to a (regular lattice) Ising model with spins (sj) living 
on the links. We will see that the spin interactions are local, albeit somewhat complicated. 

The Ising-link model is analytically and computationally much simpler than either the 
Regge or simplicial gravity approaches. But is it too simple? In Section III we present 
mean field theory results on the model in three-dimensions and in Section IV we give 
corresponding Monte Carlo results. We compare to results obtained by Hamber and 
Williams for the Regge theory in 3-d. 



II. ISING MODEL FORMULATION 



In this section we discuss how to compute the Ising action corresponding to the discrete 
form of 



c(l + bSi) is no more general, as c can be absorbed into the definitions of A and k in (2). 



S = \V-^jd''x^R, (2) 

where V is the d-dimensional volume, / d'^x^/g , and R is the scalar curvature. The lattice 
is formed out of hypercubes plus face, cubic {d > 3) and hypercubic diagonals {d > 4), 
with periodic boundary conditions [2]. First we will treat two dimensions, then three. 
Four dimensions is just like three, only harder. 

Two Dimensions : Consider a triangle with link lengths /i, I2 and l^- Define /j = 1 + bsi as 
in (1). Since sf = 1 and the formula for the area of the triangle must be symmetric in the 
3 spins, the most general form for the area is 



^123 =Co + Ci(Si + S2 + S3) + C2(SlS2 + S2S3 + S3S1) + C3S1S2S3 . 



(3) 



There are only 4 possible values for the area of the triangle, corresponding to 0, 1, 2 or 
all 3 of its spins being equal to +1. Computing these 4 areas and comparing to (3) gives 
four linear equations for the Ca in terms of the parameter b. Their solution is 

32Co = 2v^(l + b^) + 3/(6) + 3g{b) 
32Ci = 4bV3 + f{b)-g{b) 
32C2 = 2V3il + b^)-fib)-gib) (4) 

32C3 = AbVs - 3/(6) + 3g{b) , 



where f{b) = |1 - b\,J{l + 36)(3 + b) and g{b) = {1 + b)^J{l - 3b){3 - b). For example, 
b = .1 gives Ci ^ .0291, C2 ~ .0039 and C3 ^ -.0008. As stated in the Introduction, it is 
seen that we must have 6 < 1/3 for the triangle areas to be real and positive. 

Since the Einstein term in (2) is a topological invariant, it is not relevant to the case of 
fixed topology being considered here. Thus, summing over all triangles and dropping the 
irrelevant constant term, the two-dimensional action is 



S = A(2Ci "^3^ + 02 Yl ^i^j +C3 Yl SiSjSk) , (5) 

i <ij> <ijk> 

where i, j and k are link labels. < ij > indicates i and j are two of three links forming 
a triangle. In this case Sj and Sj may be termed nearest-neighbor links. < ijk > means 
that i, j and k are three links which form a triangle. The C2 term is a "nearest-neighbor" 
ferromagnetic interaction. The Ci term is a magnetic field term and the C3 term is an 
additional symmetry-breaking term. 

One might hope that the continuum limit of Ising-link quantum gravity would corre- 
spond to a second-order magnetization phase transition. But with the explicit symmetry- 
breaking terms in (5), it is clear that this transition can not be from order (< s >7^ 0) to 
disorder (< s >= 0). It would have to be an < s >^ to < s >^ transition. In two 



dimensions, as expected, we found no evidence of such a transition, at least in the mean 
field theory approximation. 

Three Dimensions : We will now go on to discuss the form of the theory in three dimen- 
sions.Q In the next section we will compare numerical results in 3-d to results for the 
unconstrained Regge theory. 

Consider the labeled tetrahedron of Fig. 1. The volume Vtet is given by the formula [2] 



144^,^, = iiiiiii-m-ii+iif-mi+ii-iir-mi-ii+iif 

+ (ii-ii+mi+ii-mi-ii+ii), (6) 

where the k may be written in terms of spins s, using Eq. (1). There are only 11 
distinct possible values for the volume of the tetrahedron, corresponding to the 11a priori 
unknown constants in the most general equation for Vtet compatible with the symmetries 
of the labeled tetrahedron: 



Vtet = C'o + Ci ^ Si + C2 ^ SiSj + C3YI SiSj 
i <ij> [ij] 

+ C4 2^ SiSjSk + C's 2^ SiSjSk + Cq2_^ SiSjSk 
<ijk> {ijk) [ijk] 

1=1 i <id> [iJ] 

Here < i,j > are again two of three links that form a triangle in Fig. 1. [i,j] are the 
remaining pairs of links. < i,j,k > form a triangle, {i,j,k) share a common site and 
[i,j, k] are the remaining triplets of links. Because sf = 1, the last four terms involve 4,5 
and 6-link interactons. Evaluating (7) for each of the 11 distinct volumes results in 11 
linear equations for the Ci and C'^. As in (4) for the 2-d case, these can easily be solved 
to determine the Ci and C'^ as functions of b. The result is not particularly illuminating 
and we will not reproduce it here. It is worth noting, however, that the the volumes are 
always real and positive if we choose 6 < 3 — y/S ~ .17. In the mean field and numerical 
results described in the next two sections, b is held equal to 0.1. 

We see that after summing up the volumes of all the tetrahedrons, the volume term in (2) 
will consist of only local interactions of the spins, involving up to 6-spin interactions. 

^One reason for not dwelling on the two-dimensional theory is that there is one (though only one) 
known disagreement between 2-d Regge quantum gravity (which appears most closely related to the 2- 
link model) and continuum results [|6|, namely the critical exponents for Ising spins coupled to 2-d Regge 
gravity 0]. One of the authors (M.G., unpublished) has examined coupling Ising spins within a fixed 
invariant distance rather than those connected by a link, since the latter coupling (used in [7]) has no 
invariant meaning. But even with the invariant coupling there was still no agreement between Regge 
gravity and continuum results; despite the fact that the total area was held fixed, enough of the links got 
suHiciently small for a finite fraction of the spins to all become coupled together. As a result the free 
energy of that model is proportional to the number of degrees of freedom squared - a fatal illness. 



The second (Einstein) term in (2), —^ J cfxy/g R, takes the form [2] 

SE = kJ2hiT.^t/^ -2n) (8) 

i t/i 

where t/i denotes a tet containing the hnk i and Otji is the corresponding dihedral angle 
at link i. For the tetrahedron shown in Fig. 1, 9t/^ is given by 

cos{e^„) = ^ ml + il - il)il - {il + il - ijM + il - il)] , (9) 

where Aijk is the triangle formed by links i, j and k. The term Se can also be written in 
terms of local spin interactions, but we shall omit the details here. 



III. MEAN FIELD THEORY IN THREE DIMENSIONS 

The Ising-link model is quite accessible to mean field theory (MFT) techniques. We write 

Z = Y.exp{-pH\s]) , (10) 

{s} 

where here [3 = 1 and 



H = S = \V-^fd^x^R. (11) 



This is a functional of the spins since /j = 1 + bsi, by Eq. (1). We wish to minimize the 
free energy, 

F= <H> -S/(3 , (12) 

for the spin probability distribution function, P[s], where S is the entropy. 

The mean field approximation [8] consists of replacing the true probability distribution for 
the spins by a factorized form: 

P[s]^p{Si)p{s2)p{s3)...p{sNi) , (13) 

where A^i is the number of links. If all links were equivalent, we could write p{si) = 
H:^ =^ J2g^ p{si) = 1 and < Si > = m. But there are three different kinds of links in 
the lattice formed out of cubes with body and face diagonals: the body diagonals, the 
cube edges, and the face diagonals. Links of the same type have the same geometrical 



environment; links of different types don't. As a result we must use the more general 
distribution, 

P]{si) = ^ , (14) 

where j = 1,2 and 3 for body diagonals, cube edges and face diagonals respectively. 

A straightforward calculation allows us to determine < H > and S as functions of mi, m2 
and m^. Let 



Pv = ^2(51)^3(52)^2(53)^3(54)^1(55)^2(56) 
Pri = ^1(51)^3(52)^2(53)^3(54)^2(55)^2(56) 
Pr2 =^3(51)^2(52)^1(53)^2(54)^2(55)^3(56) (15) 

Pr3 =^2(51)^1(52)^3(53)^2(54)^3(55)^2(55) • 



We find that 



<V>=QNoJ2...J2Pvis)Vtetis) , (16) 



Sl SQ 



< 7: I d^x^/gR> = 2itNo[7 + b{mi + 3m2 + 3m3)] 

- 6NoJ2 ■■■ E ^5(5)^t/5(5)[Pv(5) + 2Pri{s) + Pr2{s) + 2Pns{s)] (17) 

and 

S = - < lnP[s] > = -No[h{mi) + 3h{m2) + 8/1(7713)] , (18) 

where Vtet and 9t/5 are given by (6) and (9) respectively, A^o is the number of lattice sites 
and h{x) = ^-^ ln{^-^) + ^^ /n(^^). s is shorthand for si, S2, ■■■, sq. Now it is a simple 
matter to numerically minimize the free energy (12) as a function of mi, m,2 and 7773. Then 

< s > is given by 

7771+37772 + 37773 

< 5 >= , (19) 

and 



^ ,9 < f d^x./q R> , „ , ^< f d^x,/q R> , ^ 

<V > <V > 

Here we follow the notation of Hamber and Williams [9]. Also the "curvature susceptibilty" 
is defined as 



_ 2 d 



< f d^x^R> . (21) 



Results . Fig. 2 shows typical MFT results for the case A = 1. There is sharp cross-over 
behavior seen in < s >, f and TZ dX k slightly negative. TZ and < J d^x^/g R > (not shown 
in Fig. 2) are related by (20) and exhibit similar behavior in this region, xr ^^^ evaluated 
using (21), by taking a numerical derivative of < j d?x^R > with increment Ak = 1. 
We see a peak in xr here, as expected from the rapid cross over behavior in 7^. Hamber 
and Williams [9] found a second order phase transition in the 3-d Regge theory exhibiting 



XR--\kc-kt' (22) 

for k < kc, with 6 = 0.80 ± 0.06, a very weak second order phase transtion. To investigate 
whether this kind of non-analyticity is seen in the MFT approximation to the Ising-link 
model, we varied AA; from 1 downward. (AA; is the increment used to take the numerical 
derivative of < J d^x^/g R >.) The behavior (22) would result in the peak of xr growing 
with Ak like {AkY~^. 6 = 0.80 would imply that the peak would grow by 58% in height 
for each factor of 10 decrease in Ak. However, for all values of A considered (up through 
A = 75), there was no increase in the peak height as AA; was decreased from 1 down to .01. 
As a result we have no evidence of (22) with 5 — 1 < in the mean field approximation to 
the 3-d Ising-link model. Nevertheless, in the next section we will present evidence Monte 
Carlo evidence for (22) with 5 — 1 < at large values of A. 

For A 7^ 1 the situation looks similar to that shown in Fig. 2. There is a finite peak in 
Xr at k slightly negative and a first order phase transition at large positive k. The height 
of the peak in xr and the size of the discontinuity at the first order phase transition vary 
appreciably with A. Fig. 3 shows a dashed curve in A - A; space where there is a peak 
in Xr arid a solid curve of first order phase transitions. For k to the right of the first 
order phase transition, 7^ rapidly approaches zero from below as the system approaches 
a state with mi = m^ = —m2 = 1. As A; -^ — oo, the system approaches a state with 
mi = m2 = —m^ = 1. 



IV. NUMERICAL RESULTS IN THREE DIMENSIONS 



The 3-d Ising-link model was also analyzed by the Monte Carlo method. The discrete form 
of (2) was used, with toroidal topology and without higher derivative terms, as discussed 
in the previous section. As in that section, b was chosen to be 0.1 thoughtout. (See Eq. 
(I)-) 

Since there are only two possible lengths per link, there are at most 2^ = 64 possible 
configurations for a particular tetrahedron. (Actually there are significantly fewer due to 
symmetry.) Because all terms in the action are determined solely by the link lengths in 
the lattice, the limited number of distinct tetrahedrons allows many of the calculations to 
be performed only once at program entry and stored for later use in the form of "look-up 
tables". These tables are accessed during the Monte Carlo updating. As a result, the 
Ising-link model proved to be quite computationally efficient; run times were reduced by 
as much as a factor of ten over continuous-link (Regge gravity) simulations. 



Except for the reliance on look-up tables, the simulations were carried out in the usual way. 
An initial random configuration of link lengths is chosen, generating a particular initial 
geometry. A link update consists of choosing a particular link in the lattice, calculating 
the change in the action if the link takes on its other possible value, and accepting the new 
link value with probability proportional to the exponential of the negative change in the 
action (heat bath). Link updates are performed for each link in the lattice; this constitutes 
one sweep. The quantities of interest are calculated after each sweep of the lattice, and 
the values for each new geometry are binned for statistical analysis. Runs of up to 100k 
sweeps on the 4'^ and 8^ lattices, and 80k sweeps on the 16^ lattice were performed for 
various values of A and k. 

The two physical quantities of greatest interest were TZ and xrj defined in Eqs. (20) and 
(21). Hamber and Williams found that in the Regge theory, the curvature susceptibility 
diverges at points where 7?. vanishes [9]. Thus, a portion of the curve 7^ = was first 
identified (Fig. 4), and the behavior of the model was studied along that curve. But 
peaks in xr did not appear along the 71 = curve, but rather were located at values of 
k that correspond to local inflection points in TZ. This behavior is consistent with Eq. 
(21), which relates xr to the first derivative of 7?. with respect to k. Also plotted in Fig. 
4 is a dashed curve of peaks in xr- The case A = 1 was studied for values of k close to 
that curve, but no growth of the peak with increasing system size was observed, indicating 
the absence of a second-order transition in this region of parameter space. At A = 75, 
however, we did find growth in the xr peak with increasing system size indicative of a 
second order phase transition (see below). We expect that the dashed curve becomes a 
line of second order phase transitions somewhere between A = 1 and A = 75. 

Fig. 4 may be compared with Fig. 3 determined by MET. The location of the peak in xr 
is the same in both plots to within Monte Carlo statistical errors. However no second order 
phase transition occured in the MET approximation for any value of A. The agreement 
for large k between MET and Monte Carlo is poorer. As discussed in the previous section, 
MET exhibited a first order phase transition along the solid curve of Fig. 3, and for k 
to the right of that curve, TZ asymptotically approached from below. But no first order 
phase transition was found in Monte Carlo, and TZ went from negative to positive values 
at the location of the solid curve in Fig. 4. 

The comparison with MET is also seen by comparing Monte Carlo A = 1 data displayed 
in Fig. 5 with corresponding MET data shown in Fig. 2. Note the difference in the scales 
for < s > and xr- For k negative the agreement is quite good; in fact the Monte Carlo 
results agree completely with MET at A; — ;► — cxd; both indicate that the system freezes into 
a state with the body diagonals and cube edges long (s = 1) and the face diagonals short 
(s = —1). The main negative k disagreement occurs at the peak in xr which is much 
lower in the MET approximation than even on a 4^ lattice. For k positive the agreement 
is much poorer. For large k in MET, TZ never goes positive, a spurious first order phase 
transition is predicted {k ^ 48) and < s > is off by about a factor of 3. 

We now return to the evidence for a second order phase transition at large A. The desired 
signature for critical behavior would be 



n ^ Uq + A\kc - k\^ (23) 

and 



XR^B\K-kt^ (24) 

for k ~ kci where k^ is the critical point for fixed A and 5 is the critical exponent charac- 
teristic of the transition [9]. 

Previous work on the full Regge theory in three dimensions at A = 1 determined that 
TZq ~ and 5 = .80 ± .06 [9]. In that theory k had to approach kc from below; the theory 
was sick for k > kc- Here we found that the Ising-link model has TZq 7^ at the peaks 
in xr- The curvature susceptibility, though, does show behavior expected of a second- 
order phase transition. Near criticality there is an observed narrowing in the curvature 
susceptibility and an increase in peak height with increasing system size. This is readily 
seen in Fig. 6. Following Hamber and Williams [9], the finite-size scaling relation for the 
peak of the curvature susceptibility is 

a 
Hxr) --c+-lnL , (25) 

V 

where L is the system length and ajv = d{l — S)/{1 + 6), with d = 3. Using this relation, 
the critical exponent was determined from the curvature susceptibility data (Fig. 6) to be 
6 = .55 ± .02 for the Ising-link model with A = 75. 

We conclude that at least at this one value of A, the Ising-link model does appear to 
undergo a second order phase transition of the form considered by Hamber and Williams. 
However the two values of S are statistically incompatible, indicating two distinct univer- 
sality classes for the two models. 



V. DISCUSSION AND SUGGESTIONS FOR FURTHER WORK 



In three dimensions we have uncovered critical behavior in the Ising-link model for A = 75 
but not for A = 1. The region between should be carefully investigated. At A = 75 the 
critical behavior takes the form (23) and (24) as found by Hamber and Williams for the full 
Regge theory in 3-d. But T^o was in their model and not in ours. Their model was sick 
for k > kc, ours was not. At A = 75 our critical exponent 6 is statistically different from 
that seen by Hamber and Williams at A = 1 (where we found no phase transition). Does 
this difference persist at other values of A or do both models exhibit universality in A? Is 
there universality in b for the Ising-link model? Clearly more work in three dimensions is 
needed. 

The Ising-link model can and should be investigated in four dimensions as well. The mean 
field approximation is somewhat more difficult but still quite feasible, and it should be 



more accurate in four dimensions than in three. Other analytic methods should also be 
considered. 4-d Monte Carlo computations can still be performed using look-up tables 
and as a result should be many times faster than for the full Regge theory. If different 
universal behavior between the two theories persists in four dimensions, more work would 
be needed to determine which, if either theory is relevant to the universe we live in. 



ACKNOWLEDGEMENTS 



M.G. is very grateful to H. Hamber for many helpful discussions. This work was supported 
in part by the National Science Foundation through Grant No. NSF-PHY-9007497, and 
California State University. 



Note added: While writing up this work we received two related preprints by W. Bierl, 
H. Markum, and J. Riedler [10]. They independently define the Ising-link model and 
investigate its behavior in two dimensions. 
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Figure Captions 

Fig. 1. A labeled tetrahedron. 

Fig. 2. A plot of the indicated quantities versus /c for A = 1 as calculated in MFT. v is 
the average volume per site, <V > /Nq. 

Fig. 3. The 'phase diagram' of the 3-d Ising-link model in the MFT approximation. The 
dashed curve shows the location of the peak in xr and the solid curve is a first order phase 
transition. 

Fig. 4. The 'phase diagram' of the 3-d Ising-link model as determined by Monte Carlo 
simulations. The dashed curve shows the location of the peak in xr which is found to 
scale with system size at large A. The solid curve is 7?. = 0. Lattice sizes up to 16^ were 
used, and error bars are of order the size of the data points. 

Fig. 5. Monte Carlo data for A = 1 on a 4^ lattice. Data was taken every Ak = 1, and 
error bars are 0.013 or smaller. 

Fig. 6. The peak in xr for lattices of length 4, 8 and 16. The largest statistical errors of 
the data points are respectively, 0.3, 1.3 and 3.3. 



10 



